Improved nuclear matter calculations from chiral low-momentum interactions 
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We present new nuclear matter calculations based on low-momentum interactions derived from 
chiral effective field theory potentials. The current calculations use an improved treatment of 
the three-nucleon force (3NF) contribution that includes a corrected combinatorial factor beyond 
Hartree-Fock that was omitted in previous nuclear matter calculations. We find realistic saturation 
properties using parameters fit only to few-body data, but with larger uncertainty estimates from 
cutoff dependence and the 3NF parametrization than in previous calculations. 
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Major advances in nuclear structure theory over the 
last decade have been made by expanding the reach of 
few-body calculations that use microscopic interactions 
between nucleons. This progress has unambiguously es- 
tablished the quantitative role of three-nucleon forces 
(3NF) for the ground state and spectra of light nuclei 
(A ^ 12) [HQ- Pioneering extensions to larger nuclei 
reveal new facets of the 3NF, such as its role in deter- 
mining the location of the neutron dripline @,JH] and in 
elucidating the doubly-magic nature of 48 Ca [5j. Push- 
ing these successes to still heavier nuclei, which includes 
most of the table of nuclides, is a fundamental challenge 
for low-energy nuclear physics. 

The historical route to heavy nuclei is through infi- 
nite nuclear matter, a theoretical uniform limit that first 
turns off the Coulomb interaction, which otherwise drives 
heavier stable nuclei toward an imbalance of neutrons 
over protons and eventually instability. However, pre- 
dicting nuclear matter based on microscopic internucleon 
forces has proved to be an elusive target. In particular, 
few-body fits have not sufficiently constrained 3NF con- 
tributions around saturation density such that nuclear 
matter calculations are predictive. Nuclear matter satu- 
ration is very delicate, with the binding energy resulting 
from cancellations of much larger potential and kinetic 
energy contributions. When a quantitative reproduc- 
tion of empirical saturation properties has been obtained, 
it was imposed by hand through adjusting short-range 
three-body forces (see, for example, Refs. [y,0]). 

The lack of progress toward controlled nuclear mat- 
ter calculations has long been hindered by the difficulty 
of the nuclear many-body problem when conventional 
nuclear potentials are used. The present calculations 
continue an alternative approach to nuclear matter us- 
ing soft Hamiltonians derived from interactions fit only 
to few-body (A ^ 4) data. We find realistic satura- 



tion properties within our theoretical uncertainty bounds 
without adjustment of parameters. This progress fol- 
lows by applying several recent developments: systematic 
starting Hamiltonians based on chiral effective field the- 
ory (EFT) y|, renormalization group (RG) methods Q 
to soften the short-range repulsion and short-range ten- 
sor components of the initial chiral interactions so that 
conver genc e of many-body calculations is vastly accel- 
erated 13" 12' an d a new 3NF fitting procedure to the 
4 He radius rather than the binding energy Q. (Alterna- 
tive expansions using chiral interactions are described in 
Refs. [s, 13, HI)- The calculations here also employ an 



improved treatment of the 3NF contribution in many- 
body perturbation theory compared to Refs. [13, [la ], 
which includes the full treatment of 3NF double-exchange 
diagrams and corrected 3NF combinatorial factors be- 
yond Hartree-Fock. Note that previous calculations of 
neutron matter [3, E3] and finite nuclei [H, [j| arc not 
affected. 

Our results arc summarized in Fig. [TJ which shows the 
energy per particle of symmetric matter as a function of 
Fermi momentum fcp, or the density p — 2kp/(3iT 2 ). A 
grey square representing the empirical saturation point 
is shown in each of the nuclear matter figures. Its 
boundaries reflect the ranges of nuclear matter satura- 
tion properties predicted by phenomenological Skyrme 
energy functionals that most accurately reproduce prop- 
erties of finite nuclei [l^|. Although this determination 
cannot be completely model independent, the value is 
generally accepted for benchmarking infinite matter. In 
the future, calculations of the properties of finite nuclei 
will allow one to compare directly to experimental data. 

The calculations of Fig.[T]start from the N 3 LO nucleon- 
nucleon (NN) potential (EM 500 MeV) of Ref . >9]. This 
NN potential is RG-evolved to low-momentum interac- 



tions Mow*: with a smooth 
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FIG. 1: (Color online) Nuclear matter energy per particle versus Fermi momentum kp at the Hartree-Fock level (left) and 
including second-order (middle) and third-order particle-particle/hole-hole contributions (right), based on evolved N 3 LO NN 
potentials and 3NF fit to E3 H and ?"4 He . Theoretical uncertainties are estimated by the NN (lines)/3N (band) cutoff variations. 



each cutoff A, two couplings that determine the shorter- 
range parts of the N 2 LO 3NF [2l| are fit to the 3 H binding 
energy and the 4 He matter radius using exact Faddeev 
and Faddeev- Yakubovsky methods as in Ref. 22]. Our 
3NF central fit values are given in Table Q] we estimate 
that Cjj has an uncertainty of approximately 0.4 due to 
the uncertainties of the charge radius in 4 He. We use a 
3NF regulator of the form exp[-((p 2 + 3/4q 2 )/A 2 NF )"-p] 



with 



= 4, where the 3N cutoff A s 



is allowed 



to vary independently of the NN cutoff, which probes 
the sensitivity to short-range three-body physics. The 
shaded regions in Fig. [I] show the range of results for 
2.0 fin" 1 < A 3NF < 2.5 fin" 1 at fixed A = 2.0 fm" 1 . 

Nuclear matter is calculated in three approxima- 
tions: Hartree-Fock (left), Hartree-Fock plus second- 
order contributions (middle), and additionally summing 
third-order particle-particle and hole-hole contributions 
(right). The technical details regarding the treatment of 
the 3NF and the many-body calculation are as for neu- 
tron matter in Ref. [16J. We first construct a density- 
dependent two-body interaction from the 3NF by sum- 
ming one particle over occupied states in the Fermi sea 
(see also Ref. This conversion simplifies the many- 

body calculation significantly and allows the inclusion 
of all 3NF double-exchange terms beyond Hartree-Fock, 
which were only approximated in Refs. fiol [IBj ]. Fur- 
thermore, we have corrected the combinatorial factors 
at the normal-ordered two-body level of the 3NF from 
1/6 to 1/2 in diagrams beyond Hartree-Fock used in 
these references (see Refs. [1, HH for detailed discus- 
sions of these factors, which are correctly included in 
Refs. & ! E To our knowledge, previous calcu- 

lations in the literature of nuclear matter using normal- 
ordered 3NF contributions need the same correction. 

The dashed lines in the left panel of Fig. [T] (for A = 



1.8 MeV and 2.8 MeV) show the exact Hartree-Fock en- 
ergy in comparison with the results obtained using the 
effective two-body interaction (solid lines). The excel- 
lent agreement supports the use of this density-dependent 
two-body approximation for symmetric nuclear matter. 
For the results beyond the Hartree-Fock level we use full 
momentum-dependent single-particle Hartree-Fock prop- 
agators. We have checked that the energies obtained us- 
ing a self-consistent second-order spectrum overlap with 
the band of curves in Fig. [1] 

The Hartree-Fock results show that nuclear matter is 
bound even at this simplest level. A calculation without 
approximations should be independent of the cutoffs, so 
the spread in Fig. Q] sets the scale for omitted many-body 
contributions. The second-order results show a signif- 
icant narrowing of this spread over a large density re- 
gion. It is encouraging that our results agree with the 
empirical saturation point within the uncertainty in the 
many-body calculation and omitted higher-order many- 
body forces implied by the cutoff variation (the greater 
spread compared to Ref. [lij is mostly attributable to the 
corrected combinatorial factor). We stress that the cut- 
off dependence of order 3 MeV around saturation density 
is small compared to the total size of the kinetic energy 
(fa 23 MeV) and potential energy (« -38 MeV) at this 
density. Moreover, the cutoff dependence is smaller at 
fcp ~ 1.1 fin , which resembles more the typical densi- 
ties in medium-mass to heavy nuclei (p = 0.11 fm" 3 ). For 
all cases in the right panel of Fig. Q] the compressibility 
K = 175-210 MeV is in the empirical range. 

The inclusion of third-order contributions gives only 
small changes from second order except at the lowest 
densities shown. This is consistent with nuclear mat- 
ter being perturbative for low-momentum interactions, 
at least in the particle-particle channel [l(| ■ The differ- 
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ence at small densities is not surprising: the presence of 
a two-body bound state necessitates a nonperturbative 
summation in the dilute limit. We note that below satu- 
ration density, the ground state is not a uniform system, 
but breaks into clusters (see, for example, Ref. [24j). 

In chiral EFT without explicit Deltas, 3N interactions 
start at N 2 LO [2l[ and their contributions are given dia- 
grammatically by 

X 

Ci,C 3 ,C 4 C D C E 

We assume that the Cj coefficients of the long-range 
two-pion-exchange part are not modified by the RG. At 
present, we rely on the N 2 LO 3NF as a truncated "ba- 
sis" for low-momentum 3N interactions and determine 
the cd and eg couplings by a fit to data for all cut- 
offs (22j . In the future, fully evolved three- and four-body 
forces in momentum space starting from chiral EFT will 
be available (see Ref. [25| for an application of evolved 
3NF in a harmonic-oscillator basis). The fit values of 
Table U are natural and the predicted 4 He binding en- 
ergies are very reasonable. We have also verified that 
the resulting 3NF becomes perturbative in A = 3,4 (see 
also Refs. [lfl, O [13| ) , i.e., the calculated individual 3NF 
contributions are largely unchanged if evaluated for wave- 
functions using NN forces only. 

The evolution of the cutoff A to smaller values is ac- 
companied by a shift of physics. In particular, effects 
due to iterated tensor interactions are replaced by three- 
body contributions. The role of the 3NF for saturation is 
demonstrated in Fig.[5J The two pairs of curves show the 
difference between the nuclear matter results for NN-only 
and NN plus 3N interactions. It is evident that satura- 
tion is driven by the 3NF [lilQl]. Even for A = 2.8 far 1 , 



which is similar to the lower cutoffs in chiral EFT poten- 
tials, saturation is at too high density without the 3NF. 
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2.0/2.0 (EGM Ci's) 


-4.381 


-1.126 


-4.828 


-1.152 


2.0/2.0 (PWA Ci's) 


-2.632 


-0.677 


-3.007 


-0.686 



TABLE I: Results for the cd and ce couplings fit to i?3 H = 
— 8.482 MeV and to the point charge radius r4 Hc = 1.464 fm 
(based on Ref. 0) for the NN/3N cutoffs and different 
EM /EGM /PWA Ci values used. For VWk (SRG) interac- 
tions, the 3NF fits lead to Ea Hc = -28.22 ... - 28.45 MeV 
(-28.53... - 28.71 MeV). 
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FIG. 2: (Color online) Nuclear matter energy of Fig. [TJ at the 
third-order level compared to NN-only results for two repre- 
sentative NN cutoffs and a fixed 3N cutoff. 



Nevertheless, as in previous results 



13, El, the 3N con- 
tributions and the cd,ce fits are natural, and the same 
is expected for the ratio of three- to four-body contribu- 
tions. 

The smooth RG evolution used in the results so far is 
not the only choice for low-momentum interactions. A 
recent development is the use of flow equations to evolve 
Hamiltonians, which we refer to as the Similarity Renor- 
malization Group (SRG) [27M29I ] . The flow parameter 
A, which has dimensions of a momentum, is a measure 
of the degree of decoupling in momentum space. Few- 
body results for roughly the same value of SRG A and 
smooth Mowfe A have been remarkably similar (see, for 



example, Ref. [ill]). With either RG method, we can also 
change the initial interaction. The results presented so 
far all start from a chiral EFT potential at a single or- 
der with one choice of EFT regulator implementation and 
values. There are several alternatives available 
which are almost phase-shift equivalent (but the \ 2 is 
not equally good up to i?i a b « 300 MeV). Universality 
for phase-shift equivalent chiral EFT potentials as A de- 
creases was shown for smooth-cutoff V\ w k interactions in 
Ref. [tl, in the form of the collapse of different initial 
potentials to the same matrix elements in each partial 
wave channel. An analogous collapse has been found for 
N 3 LO potentials evolved by the SRG to smaller A [§]. 

Based on this universal collapse for low-momentum in- 
teraction matrix elements it is natural to expect a similar 
collapse for the energy per particle in nuclear matter. We 
consider four different chiral NN potentials: the N 3 LO 
potential by Entem and Machleidt for two different 
cutoffs 500 and 600 MeV, and the N 3 LO NN potential 
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FIG. 3: (Color online) Nuclear matter NN-only results for dif- 
ferent chiral N 3 LO potentials (EM [3 and EGM 30]). The 
upper panel shows the third-order results for V[ ow ^-evolved 
interactions at A = 2.0 fm~ (solid lines) and Brueckner- 
Hartree-Fock results for the two unevolved chiral potentials 
that provide the lowest and highest energies (dashed lines), 
EGM 600/700 MeV and EM 600 MeV. The lower panel shows 
the maximal spread of the energy results at these two cutoff 
scales A for Vi ow fc and A for SRG-evolved NN interactions. 



by Epelbaum et al. [30j (EGM) for two different cutoff 
combinations 550/600 MeV and 600/700 MeV. The re- 
sults for the energy are presented in Fig. [31 The upper 
panel shows the energies for V\ ow k NN-only interactions 
derived from different chiral NN potentials (solid lines) 
in comparison to Brueckner-Hartree-Fock (BHF, which 
means resummed particle-particle ladder) results based 
on unevolved chiral potentials (dashed lines) . For clarity, 
we only display the two extreme BHF results. As shown 
in the lower panel we find a model dependence of about 
13 MeV for the unevolved N 3 LO potentials around sat- 
uration density and approximately 2 MeV for the VWfc 
and SRG low-momentum interactions, comparable to the 
cutoff variation in Fig. [TJ The latter spread reflects the 
residual RG/SRG dependence on the initial interactions. 

By supplementing the low-momentum NN interactions 
with corresponding 3NF's we can probe the sensitiv- 
ity of the energy to uncertainties in the Cj coefficients 
(see also Refs. [3, 31, [32| )■ We consider three dif- 
ferent cases: first, we take low-momentum interactions 
evolved from the N 3 LO NN potential EM 500 MeV 
(EM a's: ci = -0.81 GeV" 1 , c 3 = -3.2GcV"\c 4 = 
5.4 GeV -1 ), second, low-momentum interactions from 
the EGM 550/600 MeV potential (EGM Cl 's: ci = 
-0.81 GeV" 1 , c 3 = -3.4 GeV" 1 , c 4 = 3.4 GeV" 1 ), and 
third, low-momentum interactions from the EM 500 MeV 
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FIG. 4: (Color online) Nuclear matter energy at the third- 
order level comparing low-momentum Vi ow fe with SRG- 
evolved chiral NN interactions for 3NF with different 
EM/PWA/EGM a values used (see text). 



potential combined with the central Cj values obtained 
from the NN partial wave analysis [33| (PWA Cj's: c\ — 
-0.76 GeV" 1 , c 3 = -4.78 GeV -1 , c 4 = 3.96 GeV" 1 ). 
The fit values for cd and eg are given in Table |TJ 

The resulting nuclear matter energies are shown in 
Fig. |4j For all three cases we find realistic saturation 
properties within the theoretical uncertainties implied 
by the cutoff dependence shown in Fig. [T] and the NN 
interaction-dependence shown in Fig. [3] The difference 
between V\ ow k and SRG results for a given set of c, is sim- 
ilar to the NN-only case (see Fig.[T]), which helps support 
the general nature of the 3NF fit. However, the present 
sensitivity study can clearly only provide a first estimate 
for the energy spread due to uncertainties of the Ci cou- 
plings. A more systematic study will require a correlation 
analysis based on a larger set of results. 

The theoretical errors of our nuclear matter results 
arise from truncations in the initial chiral EFT Hamil- 
tonian, the approximation of the 3NF, and the many- 
body approximations. Corrections to the present calcula- 
tion include higher-order many-body terms, in particular 
particle-hole corrections, and contributions from higher- 
order many-body forces and from 3NF contributions that 
cannot be expressed in terms of density-dependent two- 
body interactions. While the improvements in the cutoff 
dependence suggest that these corrections are relatively 
small, an approach such as coupled cluster theory that 
can perform a high-level rcsummation is ultimately nec- 
essary for a robust validation. 

While nuclear matter has lost to light nuclei its status 
as the first step to nuclear structure, it is still key as a 
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step to heavier nuclei and astrophysical applications like 



the structure of neutron stars [ljj . Our results can help 
with efforts to develop ab-initio density functional th eory 
(DFT) based on expanding about nuclear matter 34 1. 
This is analogous to the application of DFT in quan- 
tum chemistry and condensed matter starting with the 
uniform electron gas in local-density approximations and 
adding constrained derivative corrections. Phenomeno- 
logical energy functionals (such as Skyrme) for nuclei 
have impressive successes but lack a (quantitative) mi- 
croscopic foundation based on nuclear forces and seem to 
have reached the limits of improvement with the current 
form of functionals 
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36]. At present, the theoretical 
errors of our results, while small on the scale of the poten- 
tial energy per particle, are too large to be quantitatively 
competitive with existing functionals. The implementa- 
tion of higher-order chiral Hamiltonians and their RG 
evolution can be expected to provide more accurate and 
reliable predictions. However, there is also the possi- 
bility of fine tuning to heavy nuclei, of usin g E FT/RG 
to guide next-generation functional forms (37l l38j , and of 
benchmarking with ab-initio methods for low-momentum 
interactions. Work in these directions is in progress. 

In summary, we have presented new results for nuclear 
matter based on chiral NN and 3N interactions with RG 
evolution. The chiral EFT framework provides a sys- 
tematic improvable Hamiltonian while the softening of 
nuclear forces by RG evolution enhances the convergence 
and control of the many-body calculation. The empiri- 
cal saturation point is reproduced within our estimates 
of uncertainties despite input only from few-body data. 
Because of the fine cancellations, however, significant re- 
duction of these uncertainties will be needed before direct 
DFT calculations of nuclei are competitive. Nevertheless, 
these results are very promising for a unified description 
of all nuclei and nuclear matter. 
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